Orbital frustration at the origin of the magnetic behavior in LiNi02 
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We report on the ESR, magnetization and magnetic susceptibility measurements performed over a 
large temperature range, from 1.5 to 750 K, on high-quality stoichiometric LiNi02. We find that 
this compound displays two distinct temperature regions where its magnetic behavior is anomalous. 
With the help of a statistical model based on the Kugel'-Khomskii Hamiltonian, we show that below 
Tof — 400 K, an orbitally-frustrated state characteristic of the triangular lattice is established. This 
then gives a solution to the long-standing controversial problem of the magnetic behavior in LiNi02. 

PACS numbers; 75.40. Cx, 76.30.-v, 75.30.Et, 75.40.Mg 
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The problem of orbital degeneracy (OD) in transition- 
metal oxides has attracted considerable interest recently, 
due to its drastic effect on the magnetic, elastic, transport 
and other properties of these compounds ||l|. Among 
the most studied oxides exhibiting OD are manganites, 
vanadates, some titanates and other systems. Here we 
discuss a new interesting aspect of OD: the frustration of 
degenerate orbitals, as applied to lithium nickel oxide. 

Since LiNi02 was first studied as an ionic conductor 
the magnetic properties of this layered compound 
have been reported in numerous publications but are still 
the subject of considerable controversy. Quantum spin 
liqui d B , spin glass , frustrated antiferromagnctism 
[3M-no[7 quantum disordered state without a spin gap 
[11| and other scenarios [p2|-p0| have been claimed to 
be responsible for the magnetism of this material. Ac- 
tually, interest in this compound stems essentially from 
two facts. First, following the suggestion by Anderson 
pl| that a resonating valence bond ground state may 
exist in a S=l/2 Heisenberg triangular-lattice antiferro- 
magnet, LiNi02 was proposed by Hirakawa et al. as 
a potential physical realization of such a system. The 
second important feature, overlooked in earlier studies, 
is the unusual electronic state of Ni'^^ ions. The strong 
crystal field acting on the nickel ions brings them in the 
low-spin state {t2g, Cg] S — 1/2), with the two Cg levels 
being degenerate. Clearly, OD should play a crucial role 
in the magnetic properties of this compound ||ll| , ^ , ^ . 
However, no direct confirmation of this hypothesis exists 
yet. 

In this paper we propose a novel scenario for LiNi02 
based on an interplay between frustration on a triangu- 
lar lattice and OD of Ni^+ ions. We show experimen- 
tally that an anomalous magnetic behavior occurs in two 
different temperature regions, which we interpret as the 
consequence of two distinct energy scales present in the 



system. The first one, of the order of hundreds of Kelvin, 
is related to the antiferro-orbital coupling between Ni'^"'" 
ions and leads to a frustrated distribution of Ni'^^ or- 
bitals. The second one, of the order of tens of Kelvin, 
is due to the spin-spin interactions between electrons in 
these frustrated Ni^+ orbitals, and drives the very un- 
usual magnetic behavior at low temperatures. Our pro- 
posed scenario, which we confirm through a direct com- 
parison between experimental data and simulation re- 
sults, gives a coherent and natural explanation of this 
behavior. 

High-quality powder samples of stoichiometric LiNi02 
have been prepared following the procedure described by 
Rougier et al. [^. We have studied these samples by 
X-band {v = 9.45 GHz) electron spin resonance (ESR), 
magnetization (up to 15 T) and dc-susceptibility mea- 
surements over a wide temperature range, from 1.5 to 
750 K. 

The ESR spectra have a symmetric shape at all tem- 
peratures and show no trace of any particular feature 
that would be representative of ferrimagnetic cluster for- 
mation or of any anisotropy. This is a first but clear 
indication that the samples studied in this work are ho- 
mogeneous and of good quality. We show in Fig. 1 the T- 
dependence of the ESR linewidth, AH (half width at half 
maximum), of a typical sample of LiNi02 in the low-T 
regime. A huge increase of AH is noted when the temper- 
ature decreases below 50 K. This temperature, to which 
we shall refer below as Tm, quite naturally gives the scale 
of the magnetic interaction in LiNi02. Below T ~ 10 K, 
AH tends to a saturation value of about 0.5 T. The in- 
set of Fig. 1 shows the inverse of the spin susceptibility, 
X~^{T), deduced from the ESR spectra. As already re- 
ported earlier x{T) follows a Curie- Weiss behavior, 
X = C/{T - 0), down to about 80 K, with 6* c^i 35 ± 5 K. 
This low value of 6 confirms that the sample studied is 
so close to stoichiometry that we may consider it as pure 



1 



LiNi02 dJlltl ■ Bearing in mind that 9 represents the av- 
erage sum of exchange interactions on an atomic site, the 
positive sign of 9 imphes that the dominant interactions 
are ferromagnetic (FM). However, the Une broadening 
observed in Fig. 1 is rather indicative of strong antifer- 
romagnetic (AF) fluctuations. Line narrowing would be 
expected if the interactions were purely ferromagnetic. 
Futhermore, between 13 K and 50 K, /S.H can be fitted 
to the empirical formula A_ff — AT~p with p = 2.65. 
Such an exponent value is known to be characteristic of 
AF materials having a strong 2D character . Another 
important point is the tendency of the linewidth to sat- 
urate below 10 K. This constitutes strong evidence that 
the AF correlations do not propagate any longer below 
this temperature. 

Fig. 2 shows the magnetic moment, M , as a function 
of magnetic field, H. The curves deviate from linearity 
below Tm — 50 K, the same temperature below which 
the ESR linewidth was seen to markedly increase. At 
low temperatures, saturation is almost reached in fields 
of 15 T and, furthermore, below about 10 K, tempera- 
ture has hardly any effect on magnetization. Although 
the initial susceptibility increases with decreasing tem- 
perature in a manner reminiscent of a FM behavior, the 
M{H) curve at the lowest temperature remains below 
what would be given by a Brillouin function. This con- 
firms the existence of rather strong AF interactions along 
with ferromagnetism. 

The second anomalous aspect of LiNi02 is found at 
high temperatures and plays a crucial role in understand- 
ing the magnetic properties of this compound. Fig. 3a 
shows the results of the dc-susceptibility measurements, 
plotted in the form (T -f- lO)x(r) between 300 K and 
750 K. The curve presents a weak anomaly at around 
400 K. At higher temperatures, 9, which was found to be 
around 35 K from the low-T measurements, decreases to 
finally settle at the negative value, 9 — —10 ± 5K, for 
T ^ 600 K. The inset of Fig. 3a details the variation of 
9{T), as deduced from x{T)- Such a drastic change in 9 is 
quite unusual and gives strong evidence for a concomitant 
modification of the Ni'^^ orbital distribution in the Ni02 
layers. Additional support for this argument comes from 
the ESR data. Fig. 3b presents the T-dependence of AH, 
between 300 K and 625 K. The divergence of AH above 
room temperature is quite in line with what is expected 
when, due to changes in the Ni^+ orbital distribution, 
strong fluctuations occur. 

In order to explain this extremely rich and unusual 
magnetic behavior, we use the Kugel'-Khomskii (KK) 
spin-orbital model |jl[. Derived from the two-band Hub- 
bard Hamiltonian in the limit of a large Coulomb re- 
pulsion, this model has proven its efficiency to describe 
the interplay between magnetic and orbital ordering . 
Here we rewrite the KK Hamiltonian in the following 
form : 



H= {-J°'ri'r3 - S.Sj [Jso5{TiTj + Jdo5{TiTj + 1)]} 

<i,j> 

(1) 

where the sum runs over nearest neighbor (NN) atoms 
and 5{TiTj ± 1) is equal to 1 if TiTj ±1 = and to 
otherwise. As usual, is the spin on atomic site i while 
the pseudo-spin, t^, is introduced to describe the orbital 
occupancy. The values Ti = +1 and Ti = —1 respectively 
correspond to the two degenerate Cg orbitals of the Ni'^+ 
ion, — y^) and |3z^ — r^). In the first term, J° is the 
orbital coupling constant. J so and J do are the magnetic 
exchange constants between NN, respectively with the 
same and with a different orbital occupancy. The spin- 
orbital model is here taken in its simplest form. First, 
we consider as classical Ising spins. As a consequence, 
possible quantum effects at low temperatures cannot be 
described. Second, we do not distinguish between differ- 
ent relative orientations of the same orbital in the same 
site, as it could be done using a Potts model. Finally, we 
assume that the orbital coupling between two ions with 
the same orbital Tt = Tj = 1 is the same as the one be- 
tween two ions with — tj = — 1 . These are rather crude 
approximations but we shall see that this minimal model 
is apparently able to capture the essential physics of the 
system, at least for temperatures greater than ~ 5K. 

Recently, van den Bossche, Zhang and Mila |Q have 
suggested that the SU{A) symmetry {Jso — Jdo = J°) 
could be at the origin of the magnetic properties of 
LiNi02. Our high-temperature experimental data pre- 
sented above suggest that \J°\ » \Jso\, \Jdo\- There- 
fore, we are here testing the relevance of the alternative 
SU{2) * SU{2) symmetry in which the magnetic and or- 
bital interactions have different orders of magnitude. 

Concerning the sign of the orbital coupling, a ferro- 
orbital constant {J° > 0) would lead to a structural 
ferrodistortive phase transition. This possibility is com- 
pletely ruled out by the available X-ray data so that 
we have to consider an antiferro-orbital coupling con- 
stant (J° < 0). Since the Ni'^^ ions form a triangular 
arrangement, we then expect orbital frustration to oc- 
cur below a temperature To/ (To/ = (3/2) \ J°\ /fc^, us- 
ing the mean-field approximation) in a manner similar to 
the well-known AF-spin-frustration phenomenon. From 
the experimental results shown in Fig. 3, To/ — 400 K, 
which gives \ J°\/kB ^ 270 K. The orbitally-frustrated 
state may be described as a Wannier state (WS) [ p8[ . 
A particular property of a WS is that three NN cannot 
have the same orbital occupancy. This is illustrated in 
Fig. 4 (a and b) where we represent typical orbital distri- 
butions above and below Tof. It is evident that at high 
temperatures the probability to find a link between NN 
with a different orbital occupancy is 1/2 and that this 
probability is 2/3 in the WS. The mean-field values of 9 
are then given by: 

dht = (3./<io + SJso)/kB and 6l„,s = (4Jdo + 2Jso)/kB 
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where Bht and 9ws denote the values in the high- 
temperature state and in the WS, respectively. The as- 
sumption of AF-orbital coupling, together with the ex- 
perimentally inferred 9ht — — 10 ± 5 K and 9ws = 35 ± 5 
K, allows us to determine the two magnetic coupling con- 
stants, Jdo and Jso- We find a FM exchange between NN 
with different orbitals {Jdo/ks — 18±6 K) and an AF one 
between NN with similar orbitals (Jso/ks — — 20±6 K). It 
is interesting to note that this simple treatment of Hamil- 
tonian (1) is sufficient for the qualitative understanding 
of the magnetic behavior of LiNi02 in a large tempera- 
ture interval (T > T^). It can also be shown that these 
values of Jdo and Jso constitute the only combination of 
magnetic coupling constants that ensures the absence of 
a trivial long-range order, as indicated by the saturation 
of AH{T) at low temperatures (Fig. 1). 

To test the ability of the model to describe the low- 
temperature magnetic properties, we have performed 
Monte-Carlo (MC) simulations. The simulated system 
contains iV,; = 36x 36 coplanar ions and the classical 
Metropolis algorithm is used. At each step both 
a spin and a pseudo-spin are chosen randomly and at- 
tempts to flip them are performed independently. A MC 
step consists of Ni such attempts. After a first run of 
2x 10** MC steps to ensure equilibrium, the magnetiza- 
tion and susceptibility are computed from uncorrelated 
configurations recorded during a second run of 5x 10^ 
MC steps. The best fit to the experimental results is 
obtained for Jdo/ks = 21.0 K and Jso/ks = — 22.5K. 
Considering the simplicity of the model, the slight dif- 
ference between these values and those deduced experi- 
mentally is not surprising. The good overall agreement 
found with the experimental data is rather satisfactory 
(Figs. 1 and 2) and allows us to interpret the magnetic 
behavior at intermediate and low temperatures. In the 
range r,„ < T < Tof, the magnetic state is a param- 
agnetic one with 6 > 0. There are two thirds of NN 
pairs with different orbitals and hence a majority of FM 
exchanges {Jdo > 0) so that the system behaves as if ap- 
proaching a FM transition. Since \Jso\ > \Jdo\, at lower 
temperatures the links between the NN having identical 
orbitals are preferably satisfied (see Fig. 4c), inducing 
AF-like fluctuations because Jso < 0. This situation fi- 
nally leads to a magnetically frustrated state in LiNi02 
at low temperatures. 

To conclude, an important outcome of our experimen- 
tal results is that there are two distinct energy scales 
characteristic of the magnetism of LiNi02, which respec- 
tively correspond to the antiferro-orbital coupling be- 
tween Ni'^+ ions in the Ni02 layers and to their spin-spin 
interactions. The adjunction of two facts, orbital degen- 
eracy of the Ni'^"'" ions, and their triangular arrangement, 
leads to the build-up of a Wannier orbitally-frustrated 
state below Tof ~ 400 K. This uncommon orbital state 
is at the origin of the observed low-T complex magnetic 
behavior of LiNi02. 
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FIG. 1. Temperature dependence of the ESR linewidth, 
Aif. The straight hnc is a fit to the law ATi" = AT'P with 
p = 2.65. The inset shows the temperature dependence of 
the inverse of the spin susceptibility, as deduced from 
the ESR spectra. Circles and continuous line respectively 
represent experimental data and numerical simulations. 
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FIG. 2. Magnetic moment, M, as a function of magnetic 
field, H, up to 15 T. Circles: experimental data between 
300 K and 1.54 K; continuous curves: simulation results; dot- 
ted curve: M(H) calculated at 1.5 K using the Brillouin func- 
tion. 
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FIG. 3. a) Results of dc-susceptibility measurements, plot- 
ted in the form (T-|-10)x(T).NotethcpIatcauforT>600K. 
The inset shows the variation of 6 with temperature; the line 
is a guide to the eye. b) Temperature dependence of AH 
between 300 K and 625 K. 
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FIG. 4. Simulated distribution of orbitals at high (a) and 
low (b) temperatures. Filled and open circles distinguish be- 
tween the two types of orbitals. c) Satisfied FM links in the 
low-T Wannier state represented in (b). 
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